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1 Introduction 



Bohigas, Giannoni and Schmit JT| have conjectured that the energy levels of quantized 
chaotic Hamiltonian systems are distributed like eigenvalues of random matrices. The 
choice of the correct matrix ensemble depends on the behaviour of the classical dynamics 
under time-reversal: the unitary ensemble (CUE) for systems without time-reversal in- 
variance, the orthogonal or symplectic ensemble for systems with time-reversal invariance 
with integer or half-integer spin, respectively Although there is much numerical evidence 
to support the conjecture for systems following CUE and COE statistics, relatively few 
studies are concerned with systems with half-integer spin 0, [| [|. || . 

We shall here propose a simple model - a quantized map of the torus coupled to a spin 
precessing in a magnetic field - which is easily accessible both from the numerical and the 
analytical point of view: Suppose our particle's initial position is go with momentum p . 
With (p , Qo) fixed, we allow the spin to precess for one time unit in the magnetic field 
B(qo); we then apply the Anosov map 

0:T 2 ^T 2 , (P°)^(PA modi, (1.1) 



keeping the spin fixed; thereafter the spin is allowed to precess again for one time unit in 
B(qi), and so on. 

The idea here is that, due to the chaotic dynamics of <f>, the values of the magnetic 
field at the iterated positions q , qi, q 2 , ■ ■ ■ are essentially uncorrelated and the spin pre- 
cession therefore behaves like a random walk in SU(2). The ergodic properties of such 
constructions, so-called skew products, are well understood || [7|, §, |9f . 

The quantization of our model in terms of unitary 2N x 2N matrices U is straightfor- 
ward (section ||). In section |^ we discuss in detail those anti- unitary operators T which 
correspond to time-reversal symmetries in that T -1 UT = U _1 . If T is the only such sym- 
metry present, and if in addition T 2 = — 1 2 jv, one expects the spacing distribution for the 
eigenvalues of U to converge in the semiclassical limit N — > oo to that of the symplectic 
ensemble CSE. 

Indeed, for generic Anosov maps (p we numerically observe a transition from COE to 
CSE statistics when the magnetic field is switched on (section ^). It is, however, remarkable 



that even unperturbed cat maps, whose spectrum is highly degenerate flO|, |Tl|, exhibit CSE 
statistics for non-zero magnetic fields. 

The remainder of our paper is devoted to the semiclassical analysis of spectral correla- 
tion densities. In section || we derive a semiclassical trace formula which expresses sums 
over eigenvalues of U in terms of sums over fixed points of the classical Anosov map (p. We 
observe that the fixed points remain unchanged when switching on the magnetic field but 
that additional weight factors, taking care of the spin contribution, appear. In the case of 
flows an analogous observation has been made in |T^, 13|. In section the trace formula is 



employed to investigate two-point correlations of the eigenvalues of U, with focus on the 
behaviour of the spectral form factor for small arguments. In particular, we show that the 
diagonal approximation |TJ], [T5|] reproduces the random matrix prediction if the magnetic 
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field is sufficiently generic. We essentially follow the line of arguments in fTE |, where the 
form factor for continuous-time dynamics is discussed. However, the assumption in |16] 
that the skew product dynamics must be mixing is replaced by a simpler condition. In the 
final section [7] we discuss the interesting fact that the diagonal approximation also works 
if the spin precession does not become equidistributed in SU(2), but only takes a finite 
number of distinct values in SU(2). 



2 Maps with spin 

The Hamiltonian of a particle with spin ~ in a magnetic field B(q) is given by (in units 
where 2mc/e = 1) 

E = E l 2 -h(T- B{q), (2.1) 

where H is a scalar Hamiltonian, and the second term describes the interaction between 
the magnetic moment of spin and the magnetic field B(q). Here I2 denotes the 2x2 unit 
matrix, and er is the vector of Pauli matrices 

*«=(J J), *y=(i Q 1 )' ff *=(o -l)- (2 ' 2) 

We are interested in modeling free spin precession interrupted only by periodic delta- 
shaped kicks at integer times t = n G Z. Such a system is represented by the time 
dependent Hamiltonian 

K(t) = ^2K ± 2 5(t-n)-htr-B(q), (2.3) 

cf. j|, H]. The Floquet operator which propagates the system, e.g., from time t = — ~ to 
t = ~, is then given by 

U = U s U t U s , with U s = e^ B{q) and U t = e -^ H °l 2 . (2.4) 

In this paper, the translational dynamics U t will be represented by a quantized map of 
the torus T 2 , which is given by a unitary operator U^ acting on an N dimensional Hilbert 
space T~Ln — C N . The dimension of the Hilbert space 1~Ln and Planck's constant are related 
via the condition 2rchN = 1. (We refer the reader to appendix A for more details on the 
quantization of maps on the torus.) For instance, in the case of the map 



the action of Un on functions $ G H.n is given by |1 

1/2 N-l 



(-, \ 1/2 JV -1 |-„ . 

' ri' = n L 



$(Q') (2.6) 
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Q = 0, . . . , N — 1. We shall also investigate the perturbed cat map 

4> : po A o p (2.7) 

with perturbation 



p-[ p y{ v+ \ m ) -di, (2.* 



where / is a bounded periodic function and k measures the strength of perturbation. 
Since Anosov systems are structurally stable, 4> remains Anosov if k is small enough. The 
quantization of <fi is represented by [[L8| 

\J N {<f>) = Ujv(p) U N (A) XJ N (p) (2.9) 

with 

[Ujv(p)$](Q) = e* iNkS (%)<!>(Q) (2.10) 

and /(g) = 

The propagators U t (A) and U t (</>) now act on a two-spinor \1> e C 2 ® TCn simply by 

tl J VU^(A)$ 2 j' Utlw VUjv^)^' V^y' ^'^e/tiv, 

(2.11) 

and the action of U s is straightforwardly given by 

[U s tf](Q) = e^ B (#)^(Q). (2.12) 

3 Anti-unitary symmetries 

Following Bohigas, Giannoni and Schmit U , we expect that the eigenvalue statistics of the 
quantized map U converge to those of the circular random matrix ensembles, CUE, COE, 
or CSE. The choice of the correct ensemble depends on whether the quantum map possesses 
an anti-unitary symmetry T. In particular, the spectral fluctuations are expected to agree 
with those of the CSE ensemble if T 2 = — 1 2 jv- In this section we discuss the necessary 
conditions to be imposed on the classical map <f> and on the magnetic field B(q) in order 
to observe CSE statistics. 

To this end consider an Anosov map (f> invariant under time reversal, i.e. 

to0ot = </T\ where r : ( J i-» ( j . (3.1) 
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Moreover, suppose that the map <f> is also invariant under inversion, i.e. 

7r o (f) o 7r = <fi with 7T : 
Then, <fi will also be invariant under a transformation with 

f = 7T o T, (3.3) 

which we will call non-conventional time-reversal, see, e.g. []19| for a general discussion. 
For spinless motion the quantum time reversal operator K is simply complex conjugation, 




[K¥](Q)=*(Q), 
and the operator corresponding to f is given by 

[K*](Q) = *(-Q). 



(3.4) 



(3.5) 



Moreover, since r and r are symmetries of 4>, the quantum map Ut(0) will be invariant 
under K and K, i.e. 



KU t (0)K =U- 1 (0), KU t (0)K =U" 1 (0). 



(3.6) 



In general, however, the coupled map U(0) = U s U t (</>)U s is not invariant under these 
symmetries, because 



[KU s Kf] (Q) = e-3^(#>tt(Q), 
KU S K^1 (Q) = e -^-s(-#)^(g) ; 



(3.7a) 
(3.7b) 



where B(q) = (B x (q), —B y (q),B z (q)). Since we are concerned with system with spin, we 
should also invert the direction of spin when reversing time. Let us consequently define 
the modified time-reversal operators T and T by 

[Ttt](Q) = e^ 2 ^(Q) = i* y *(Q), [T *](Q) = U7„*(-Q). (3.8) 
Note that T 2 = T 2 = — I2N, as needed for CSE statistics. For these operators one obtains 



[T-^T*] (Q)=e^ B ^(Q), 
(Q)= e L 2° B (-^(Q). 



T -i U,Ttf 



(3.9a) 
(3.9b) 



From ( |3.9a| ) we easily see that T 1 U s T / U s 1 for non-zero magnetic fields, since obviously 
B(q) 7^ ~B(q). However, rel. ( p.9b|) implies that T _1 U s T = U" 1 for odd magnetic field, 
i.e. if B(-q) = -B(q). 
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In addition we have to make sure that there is no further anti-unitary symmetry, in 
particular the system must neither be invariant under K nor under K. If B y (q) vanishes 
identically then B(q) = B(q) and by (|3.7a|) the spin part is also invariant under K, 
which should be avoided. The last argument is of course independent of the choice of the 
coordinate system, and hence we require that there is no vector r ^ such that r-B(q) = 
for all q, i.e. the magnetic field B(q) must have linearly independent components. 

To summarize, in order to observe CSE spectral distributions we shall require that f 
is a symmetry of the classical map (p, that the magnetic field B(q) is odd and that its 
components are independent functions. 



4 Spectral statistics 

After having clarified under which conditions we expect to see CSE statistics we will 
now test our arguments in a numerical experiment for the quantized cat map A and its 
perturbation (f>. A simple choice of the shear in (|2.8|) is, e.g., given by 



/ ( g ) = i_ sin2v rg. (4.1) 



The map <f) will be Anosov if [18 



k < fc max = ^jJ- = 0.327... (4.2) 
V5 



The only symmetries of the cat map ( |2.5| ) are time- reversal r ( |3.1| ), inversion iz ( |3.2|) , and 



non-conventional time-reversal f ( p.3|) . The choice of the perturbation ( |4.1|) , with / being 
odd, conserves all symmetries r, f and 7r. 

Now, in order to observe CSE statistics, the spin precession must be caused by a 
magnetic field which is odd and periodic in q and whose components are independent 
functions (cf. section |3|). These constraints guarantee that the only anti- unitary symmetry 
of the quantum system is T ( ^.8| ). For our numerical calculation we chose the magnetic 
field 

(sin(27rg)\ 
sin(4vrg) , (4.3) 
sin(67rg) J 

which clearly satisfies these requirements. 



We now combine \J N (A) and Utv(</>) with the spin precession propagator U s as in (|2.4| ) 
to obtain the quantum maps 

U(A) = U s U t (A)U s and U(0) = U s U t (0)U s , (4.4) 

where 
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The matrix elements of U(A) and U(0) are easily worked out using ( |2.6| ) - ( |2.12| ). If 
N is divisible by four, there exist simple relations between the matrix elements of the 



propagators in eqs. (|2.6| ) and (SH) |20|j which lead to invariance with respect to the unitary 
operator 

[Rtf] (Q) = ie^(Q + N/2). (4.6) 

This operator obeys the relations 

[U(A),R]=0, [U(0),R] = O, [R,T] = and R 2 = -t 2N . (4.7) 

As a consequence in this case the statistics of the spectra of U(A) and U(</>) would follow 
those of the CUE ensemble, cf. [|19]. Thus, we shall require that N ^ mod 4. 

Studying the spectral statistics of quantum maps means looking at the distributions of 
the eigenangles {a\, . . . , a^jv}- In order to compare energy fluctuations in different systems, 
the spectra must be unfolded, that is the energy levels are rescaled so that the mean level 
spacing is one. (The spacing s is defined as the distance between two consecutive levels.) 
Moreover, in order to observe universal distributions the spectra must be desymmetrized, 
that is only energy levels with the same quantum numbers corresponding to all mutually 
commuting unitary symmetries must be considered. Therefore, when analyzing the spectra 
of Un{A) and U(0) in the absence of spin dynamics, we have to take into account only 
eigenphases which correspond to the same eingenvalue P = ±1 of 7r. For systems with half- 
integer spin and time- reversal invariance an additional proviso must be taken into account: 



each energy level has at least multiplicity two (Kramers' degeneracy |HJ]) because the 



square of time-reversal operators T for systems with half- integer spin being — ±2N causes 



T\& to be orthogonal to \& f22| . The statistical analysis, however, must be performed on a 
spectrum obtained by removing such degeneracy. Compare section |^ for details. 

Let us briefly describe the spectral distributions of XJn(A) and Uat(0) in the absence of 
spin dynamics. The spectra of quantum cat maps Un(A) are well known to be highly non- 
generic, and in particular do not follow any universal distribution ||23||. The reason for this 



untypical behaviour is the high number of quantum symmetries - about O(N) many |TI| |TT 



which commute with the quantum propagator Ujv(-A). Such arithmetical symmetries 



can easily be broken by slightly perturbing the classical map fllSfl , and the expected COE 
random matrix statistics are recovered. This behaviour is evident in fig. |], where the 
consecutive level spacing distribution p(s) and the integrated level spacing distribution 

p(s')ds' (4.8) 



are plotted for the quantum maps ( |2.6| ) and 

Fig. ^] shows the spectral statistics of U(A) and U(0) defined in ( |4.4| ) with the magnetic 
field given by (|4.3| ). Clearly, in both cases the spectral distributions nicely follow CSE 
statistics, as predicted. It is remarkable that the non-generic behaviour of the quantum 
cat map without spin (figs. |l(a)| and |l(b)| ) has disappeared, even though the spin coupling 
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(a) Spacing distribution of the quantum cat 
map (2.6) with P = +1. 



(b) Cumulative spacing distribution of the 
quantum cat map (2.6) with P = +1. 
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(c) S pac ing distribution of the quantum 
map (U) with k=0.32 and P = +1. 



(d) Cumulativ e sp acing distribution of the 
quantum map fl2.9|) with k=0.32 and P = +1. 



Figure 1: Local spectral statistics of the quantum cat map ( |2.6|) (above) and of its pertur- 
bation (|2.9|) (below). The eigenangles belong to states with the same eigenvalue of ir and 
their number is iV = 1021. 
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Figure 2: Local spectral statistics of the quantum cat map (above) and of its perturbation 
(below) when the spin dynamics is introduced. The number of eigenangles N is 1021. 
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is of sub-principal order in h and therefore does not affect the highly degenerate classical 
dynamics. The semiclassical arguments in the following sections will explain the CSE 
statistics for generic maps 0, but in the case of the unperturbed cat map A the situation 
is more subtle. We will briefly return to this point in section ^ when discussing piecewise 
constant magnetic fields. 



5 Semiclassical trace formula 

The purpose of this section is to derive a trace formula for the quantum maps fl4.4|) . Since 
Berry's study of the spectral rigidity [15j based on Gutzwiller's trace formula |24]], trace 
formulae are the main semiclassical tools in analyzing spectral correlations. In the case of 
flows trace formulae including spin contributions have been derived in [13] for the Pauli 
and the Dirac equation. 

In what follows <p can be any Anosov map on T 2 . The magnetic field B will be left 
generic too, except that, for simplicity, it should be a function of q only. Our goal is to 
determine Tr[U n (0)] as a sum over periodic orbits. For this purpose it is more convenient 
to consider the operator 

U(0) = U t (</>)U 2 . (5.1) 

Since U(0) is unitarily conjugate to U(0), clearly Tr[U n (0)] = Tr[U™(0)]. 

Let Op A r(/) denote the Weyl quantization of a classical observable / G C°°(T 2 ) (see 
appendix |A]). For a 2 x 2 matrix 

M=(^H ^=(M)eT 2 , (5.2) 



K 92i{z) g 2 2{z 

whose elements are smooth functions on T 2 , we then define 

Opjv(pn) Op N (g 12 ] 



which semiclassically fulfills (cf. appendix [A]) 

Op N (g k 9k-i ■ ■ -9i) ~ Op N (g k ) Op N (g k -i) ■ ■ .Op^foi), N -> oo. (5.4) 
With the choice 

g(z) = exp (icr • B(q)) G SU(2) (5.5) 

we have U 2 = Op N (g) and, thus, we may write U(</>) = U t (</>) Op N (g). We now rearrange 
U n (0) as follows: 



L>(0) = U»(0) (Uj-»(0) Op^^UT 1 ^)) (U t 2 ' n (0) O P7V (^)Ur 2 (0)) ■ ■ ■ 

• • • (U- 1 ^) Op N (g)V t (cf>)) Op N (g). (5.6) 
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Applying Egorov's theorem (appendix 

Uf-"(0) Op N (g)Vr k W ~ Pjv fo o r- k ), N^oo (5.7) 
to each factor (k = 1, . . . , n) yields 



n-l 



U B (0)~U?(0)n°P^°0*), M) 

fc=0 

where the product is time ordered, i.e. 

n-1 

J] Op^ o fc ) = Pjv fo o 0™" 1 ) Op N (g o </>™" 2 ) . . . Op N (g). (5.9) 
fe=0 

In order to simplify the notation, it is convenient to set 

g k (z)=go ( f> k - 1 (z), k = 0,...,n, (5.10) 



with the convention that go(z) = I2, and g%{z) = g(z) = exp (icr • B(q)). Now eq. (|5] 
becomes 



U w (0)~U^(0)J];Op^ fe ), iV^oo. (5.11) 



k=0 



In the case when <p = A is an unperturbed cat map, Egorov's theorem is an identity, not 

merely an asymptotic relation, and hence eq. ( |5.11|) is an identity too. Due to ( |5.4j) in the 
semiclassical limit (N — > 00) eq. ( |5.11| ) becomes 

~ U?(0) Op^), (5.12) 

where 

n 

d n (z) =Y[g k (z)e SU(2). (5.13) 

fc=0 

Let <fi, . . . , <^7v be an orthogonal basis of eigenfunctions of Uat(</>) with eigenphases Uj, 

XJ N (<l>)<p j = d u i<p j . (5.14) 
Taking the trace of the right-hand side of ( 5.12Q yields 

N 

Tr [U?(0) Op N (d n )] = VN (tvd n )^) e 1 ^. (5.15) 
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The leading order term in the semiclassical limit N — > oo of the above can be shown to be 
(cf. appendix |B|) 

V(^,0 PiV (trd n )^)e in " J ~ V] tT ^A exp [27rW (5>(gf + ™i, q f ) - m 2 q f )] . (5.16) 

Here Zf are fixed points of order n with winding numbers (m^mj) = m, i.e. 

0"(z f ) = z f + m. (5.17) 

Moreover, Ri/ = det(Ai™ — I), where M%, '■— ^-( z {) is the monodromy matrix and 

is the generating function of <p n on the covering plane. We will denote the classical action 

of the fixed point Zf by 

SV(z f ) := S ( ^{q i + m 1 ,q i ) -m 2 gf. (5.18) 
Thus, in leading semiclassical order (N — > oo) we find 

Tr [IT (</>)] ~ V ^tol exp [27riiVV(*f)] , (5-19) 

which is the main result of this section. Note that ( |5.19| ) is always a semiclassical expression, 
even if <fi = A is a cat map. Semiclassical approximations enter in ( |5.12 ) and in ( [5.16 ) 
for both the unperturbed and perturbed map. Formula ( 5.19 ) shows that the classical 
orbits are not affected by spin precession and that the spin contribution is represented by 
the factors trd n (zf) for each periodic orbit. The respective result in the case of flows was 
obtained in JTJ, |T3|. 



6 Semiclassical analysis of two-point correlations 

We will now apply the trace formula derived in the previous section to the study of spectral 
two-point correlations, essentially following the line of arguments of fl6fl , giving a new 
argument concerning equidistribution in SU(2). 

Since U is invariant under T with T 2 = — 1 2 jv, cf. section [|, it is easily seen that the 
eigenvalues must have at least multiplicity two (Kramers' degeneracy). Let us denote by 
, a 2 N ^ [0, 27r) the eigenphases of U with eigenfunctions tyj, i.e. 

VM'j e in 'vl/ r (6.1) 

labelled such that = ctk+N, for k — 1, . . . , N. We will then only need to consider 
correlations between the first N distinct eigenphases ... ,ajy. 
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In order to measure correlations on the scale of the mean level spacing, we rescale the 
spectrum by putting 

N 



2tt 



Oik, 



k 



N. 



The two-point correlation density of this sequence is defined as 



1 N 



(6.2) 



(6.3) 



k,l=l mf 



and its Fourier transform, the spectral form factor K 2 (t, N), reads 



K 2 (r,N) 



N 



R 2 (s,N)e ZmTS ds 
' 'Tr(U")| 2 -iV5 n o 



N 



N 

fc=i 



2-K\XuT 



-N5, 



(6.4) 



AN 



for r = jj, n G Z. According to the Bohigas-Giannoni-Schmit conjecture 0, one expects 
that semiclassically (N — > oo) the form factor K 2 (t, N) converges on average to the corre- 
sponding CSE form factor. That is, for any smooth and rapidly decaying test function 0, 
we expect 



lim Y, K *(%> N ) = I K^{r)<j>{r)dT, 



(6.5) 



with 25 



K. 



CSE, 



T 



~|t| — ||t| log |1 — t\ for |r| < 2 

for \r\ > 2. 



(6.6) 



This statement implies that also the pair correlation density R 2 (s,N) converges (on aver- 
age) to the corresponding CSE density. 

It seems to be extremely difficult to prove rel. ( |6.5| ) with present techniques. Here, 
we will aim at understanding the asymptotics of the form factor for small values of r, 
as iV — > oo, in the regime governed by the diagonal approximation [|14], 15| . For the 
diagonal approximation to hold, we will assume in the following that there are no systematic 
degeneracies in the classical actions of the map under consideration. This is true for generic 
perturbed cat maps 0, but not for the original cat maps A, where actions are in fact highly 
degenerate p3[ p6fl . 



By virtue of our trace formula (|5.19|) we have, for t = jt 0, N ^ oo, 



K 2 {t,N) 



— lTr(U")| 2 ~ 



AN 



1_ trd n (z f ) tr d n (z' { 



z t ,z' 



(n) 
Zf 



R^ 



exp [27riJV(V(«f3- V(*f))l 

(6.7) 
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(recall tr d n (z' { ) = trd n (z' { ) for d n G SU(2)), where the sum extends over all fixed points 
of order n. Since for S < pn(Z{) ^ S^n^z'f), cf. ( |5.18|) , the exponential in (|6.7|) shows rapid 
oscillations as N — > oo, we assume that in the combined limit 

r^O, N-+00, n = rN^oo (6.8) 

the double sum in (|6.7|) is dominated by the diagonal terms [14], p~5| , i.e. by the terms with 
S(j) n {zi) = Scj,n(z' f ). We know that 

V(*f) = V(*f), 4? = #4, trd n (z { ) = tr d n (z' { ) (6.9) 

for all points z' f along the periodic orbit (compare section Q), 

z' t = 0(z f ), . . . , ^(zf), n (* f ) = z f . (6.10) 

Thus, we find n* degenerate actions for every orbit of period n, rfi denoting the 'primitive 
period, i.e. n = krfi , k G N. Furthermore, since our system is invariant under (non- 
conventional) time-reversal T, we have 

V(*f) = S r (f(z { )), R<g = 4^ f) , tr d n (z { ) = trd n {r{z { )). (6.11) 

Neglecting self-retracing orbits in the limit n — > oo this results in an additional factor of 2 in 
the diagonal approximation. Note that if the classical map <f> has a further symmetry, e.g. 
inversion 7r ( |3.2|) , we also have S^,n(Tr(zf)) = Sf^Zf). However, since in general tr d n (z{) 
and tr d n (7r(zf)) become uncorrelated for large n, we will neglect cross correlations of these 
terms. We therefore conclude that in the combined limit (|6.8|) the form factor can be 
approximated by 

^-^^f. ,6,2, 

It is well known that the number of fixed points of Anosov maps grows, on average, 
like ~ e hn for n large, where h denotes the topological entropy. What is more, the fixed 
points become equidistributed in phase space T 2 0, [27] , in the sense that for any smooth 
test function a(z), we have for n large 



where the average (...)„ is some linear mean over an interval about n whose size grows to 
oo as n — > oo. Let us ignore the fact that d n (z) depends on n, since, as we shall justify 
below, it in fact converges on average to a constant independent of n. Hence by virtue of 
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In the semiclassical limit the quantum dynamics reduces to the skew product dynamical 



system 16 



Y : T 2 x SU(2) -> T 2 x SU(2) 

(z,g) >-> (4>(z), gi(z)g) 



(6.15) 



with gx{z) = exp (icr ■ B(q)). The rath iterate is then given by 

Y n (z,g) = (<f> n (z), d n {z)g). (6.16) 
The right-hand-side of (|6.14 ) can be viewed as a sequence of probability measures 



F{Y n {zM))) dz 



(6.17) 



if we put F(z,g) = (tig) 2 . For F bounded on T 2 x SU(2), this sequence is contained in 
a compact space of probability measures, hence we find a convergent subsequence rij with 
some limit measure /x, i.e. 

lim v n .{F) =fi{F). (6.18) 

J— »oo 

Due to the (. . . )„ average we have, for n large, 

v n (FoY r )~ Vn (F), (6.19) 
for any fixed integer r. Moreover, (with the substitution w = 4>{z)) we have 

u n (FoY)= [ (F(^ n+1 (z),d n+1 (z))) dz= [ (f (<t> n (w), d n+1 {d>- l (w)))\ dw 
[F(<t> n (w), d n (w)g 1 ((f>- 1 (w))) y^dw = v n {Fo ll ) 

(6.20) 



where j 1 is defined by 

7i : (z,g) >-> («, • 
Therefore, we have ^„(^°7i) ~ u n(F), and similarly, since 

/ F(0" +r '(2), d n+r (z)) dz= [ F(<f> n (z),d n (z)d r ((l>- r (z))) dz, 

for any fixed r G Z we obtain 

u n (Fo lr ) ~ z/ n (F), 



(6.2i; 



(6.22) 



(6.23) 
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with 7 r : (z,g) i— > (z, g d r ((/)~~ r (z))) . Thus, we conclude that the limit measure \i itself 
must be invariant, i.e. \x{F o y r ) = fi(F), for all r6Z. Clearly, since \i is invariant under 
7 r , it is also invariant under 77 1 an d i n particular we obtain the relation 

/i(Fo7 r )~/i(F) with 7r := 7r o 7r -_ 1 1 :(z^)^(2 ; , g gi(</>^ (z))) . (6.24) 

So if, for almost all z G T 2 , we find a set V of integers such that the group generated 
by {gi(4>~ u (z))} u( zv is dense in SU(2), then the action of that group is uniquely ergodic on 
SU(2) (see |28j for more details and references on equidistribution on SU(2)), and we have 
that dfj, = dzdg, where dg denotes Haar measure. Clearly it is always possible to satisfy 
the above condition for any generic choice of gx, i.e. for any generic magnetic field B(q), 
cf. section |3|. 

Since the limit /1 of every converging subsequence is unique, in fact every subsequence 
converges to \i. That is 



lim / ( F(Y n (z,t 2 )) ) dz = / / F(z,g)dzdg. (6.25) 

n->oc Jj2 \ I n Jj2 isU(2) 

In our case F(z,g) = (tig) 2 , and by the character formula [29[] 

/ / F(z,g)dzdg= [ (trg) 2 dg = l. (6.26) 

JT 2 JSV{2) </SU(2) 

Therefore, for generic magnetic fields B(q), the asymptotics of the form factor at small 
t is in the diagonal approximation given by 

K 2 (t,N)~ 1 - T) (6.27) 

which is identical to the small r asymptotics of the CSE form factor ( |6.6|) . 



7 Piecewise constant magnetic fields 

As we have seen in the previous section, the diagonal approximation works when the spin 
precession becomes equidistributed in SU(2). It is quite remarkable that this is not a 
necessary condition. Suppose for instance that the magnetic field is piecewise constant, 
such that (yf 1 (0 1 (2)), (7x((/) 2 (z)), . . . only takes values in a discrete set in SU(2). In this 
case the equidistribution theorem ( |6.25| ) holds with SU(2) replaced by T, where T is the 
finite subgroup generated by {gi(<fi r z)} r , 

lim / /f(Y»(*,1 2 ))\ d 2 4E / F{z,9)*z, (7-1) 

n-WTP \ / n \ T \ ^J^ 2 

where |T| is the order of V. Now 

^ 2 (r,iV)~iri-^(tr^) 2 . (7.2) 
1 1 ser 
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If the representation of T in SU(2) is irreducible, we have as a consequence of Schur's 
Lemma (see, e.g., |29|] for details) 



-^(trs) 2 = l, 



(7.3) 



as above. Irreducibility is guaranteed by our assumption in section |3] that there be no 
vector r ^ such that r ■ B(q) =0 for all q. Hence we have again 



(7.4) 



for r small. 

Let us complement this argument by a numerical experiment. We choose the magnetic 
field 



B(q) 



'(ti-,0,0) if0<g<l/6 

(0,vr,0) if 1/6 < g < 1/3 

(0,0, tt) if 1/3 < g < 1/2 

1(0,0,0) if g = 1/2. 



(7.5) 



The symmetry constraint B(l — q) = —B(q) determines the magnetic field in the interval 
1/2 < q < 1. We also require periodicity, i.e. B(q + m) = B(q), m e Z. This particular 
choice leads to the discrete group of Hamilton's quaternions, 

T = {±1 2 , ±ia x , ±ia y , ±\a z }. (7.6) 

The quantum maps U{A) and U{<fi) defined in ( f4.4|) can also be easily diagonalized with 
the magnetic field given by ( |7.5| ). The statistics are clearly CSE as shown in fig. [3]. This 
outcome could be expected for generic Anosov maps, since the diagonal approximation for 
the form factor agrees, at small argument, with the random matrix prediction. 

In the case of the unperturbed cat map, however, we cannot proceed as in section [5] due 
to the exponentially large number of fixed points sharing the same classical action |3(| ^3] . 
Remarkably, we still observe CSE statistics (fig. ||). Similar peculiarities are found for 
arithmetic triangles with non-conventional boundary conditions |3l|, non-arithmetic Hecke 
triangles p2| and hyperbolic tetrahedra |33| . 
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Figure 3: Local spectral statistics of the quantum cat map (above) and of its perturbation 
(below) when the spin precession is caused by the magnetic field (|7.5|) . The number of 
eigenangles N is 1021. 
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A Quantum mechanics on the torus 



We briefly review quantum mechanics of systems whose classical phase space is the torus 
T 2 . For more details see |T7], |3§. 

Because of the topology of the torus, quantum states are taken to be periodic in both 
position and momentum representation, i.e. 

i[)(q + mi) = i/i(q), ip{p + m 2 ) = V>(p), mi,m 2 eZ (A.l) 

where 



iP(p) = —= ^(q)e-^dq. (A.2) 

The periodicity of the wavefunction in both bases has two important consequences. Firstly, 
both ip(q) and ip(p) are superpositions of delta functions supported on the lattices points 
q = 2irh Q and p = 2irh P respectively, where Q, P G Z, i.e. 

^) = ££ 6 ( q -^ + m ) (A - 3) 

with *5/(Q + N) = ^(Q)- Secondly, 2irh must be an inverse integer, i.e. A" = l/2nh. It 
follows that the Hilbert space may be identified with the A^-dimensional vector space T~Cn 
with inner product 

Q mod N 

In order to quantize observables / G C°°(T 2 ) we need to introduce the translation 
operators 

ti$(Q) = $(Q + l) (A.5a) 
t 2 $(g) = $(Q)exp^Q^, (A.5b) 

which may be viewed as the exponentials of the usual differentiation and multiplication 
operators on the real line. For any m, n G Z, we have the following commutation relation 

tytl = exp ( — mn) t n 2 t™ . (A.6) 



Note that = — Ijy. The Weyl-Heisenberg operators are defined by 

rjv(n) = exp Ul n 2 \ t?t?, (A.7) 
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where n = (m, ^2)- We then have the following multiplication rule 



7T1 



T N (m)T N (n) = exp ^— (m, „)J T N (m + n), (A.8) 

where u (m, n) = m\n2 — m^Ux is the standard symplectic form. 

Let / G C°°(T 2 ) be a classical observable on T 2 whose Fourier series is given by 



f(z) = Le 2 ™ m , z = (p, q) e T 2 . (A.9) 



The Weyl quantization of / is defined as 

Op*(/)= E f™T N (rn). (A.10) 

Since Op N (f*) = Op N (f) jf , Op JV (/) is Hermitian if and only if /(z) is real. 

Semiclassically, quantum observables commute, which can be easily seen by expanding 
the exponential in ( A.8|) . More precisely, let f,g £ C°°(T 2 ), then 

Op N {f)Op N {g)^0 PN (fg), N^oo. (A.ll) 

If / and g depend only on either p or q, the above relation is an identity for each N. 
The quantization of classical matrix-valued observables is analogous to the case of scalar 
functions. Let g(z) be a 2 x 2 matrix such that g^ G C°°(T 2 ). We define its Weyl 
quantization as 

OpM = (n PNi ( gn \ n VN ( 9 i) (A. 12) 

From ( A.ll ) it follows that 

Op N (gk ■ ■ ■ 9291) ~ Op iY (^ fc )...Op 7V (^ 2 )Op 7V (^i), A ^00. (A.13) 
Clearly, we have 

OpjvG^) = Pjv (^)t. (A.14) 

If g(z) is unitary, combining ( |A.13| ) and ( |A.14| ) yields 

Op N (g)Op N ( g y ~ 1 2JV , A ^00, (A.15) 

that is Op A r(^) is only semiclassically unitary. 

We can now formulate a version of Egorov's theorem |36j for matrix valued observables. 

Let 

W)=( U " W Uk °J. (A.16) 
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where Ujv(0) is the quantum propagator of an Anosov map <f>. We have 
tt-V^Ot, (n\u (th\ - ^ u ^(^) _1 °Pjv(^i) u ^(^) U iV (0) _1 Op iV (^i 2 )Ujv(</)) 

u t w^^UtW-^u^-iop^ou^) u^^Op^u^) 

OPiV^ll ° 0p w (#i2 o 0) 

0piv(^21 O 0) Op w (>22 O 0) 

as N — > oo. For cat maps this relation is an identity. 



(A.17) 



Op N (g o <f>), 



B Trace formula for matrix elements 

The following trace formula for matrix elements of quantum observables is a well known 



generalization of Gutzwiller's trace formula |37|, |38|, p9|| . We will here present its derivation 
in the case of general Anosov maps (f> on the torus T 2 , which is particularly clean and 
simple. The formula is needed in section [|. 

Suppose <fi, . . . , ifN is an orthogonal basis of eigenfunctions of Ujv(</>) with eigenphases 



UnWw = e ia >;- (B.l) 

Claim: For any smooth classical observable / G C°°(T 2 ) and n fixed we have in the 
semiclassical limit (N — > oo) 

V(^,0 Pjv (/)^)e in ^ ~ ex P l 27TiN (V(gf + m u q f ) - m 2gf )] , (B.2) 

where (mi, m.2) = m are the winding numbers and Z{ are the fixed points of order n, i.e. 

n (z f ) = 2 f + m. (B.3) 

Furthermore, = det(.M™ f — I), where M.^ = ^(^f) denotes the monodromy matrix 
and Sc/y generates the classical map <p n . 

The left-hand-side of ( |B.2| ) is of course the trace of U^(0) Op N (f). Since 

and because / G C°°(T 2 ) implies that the coefficients / M are rapidly decreasing, we 
may without loss of generality consider only the matrix elements U^(</>)T/v(/^). For 
v = iV — 1 the functions 



w) = r ^=^ modiV ) (b.5) 

" ^0 (Q^z/modiV) 1 ; 
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form an orthonormal basis of the Hilbert space T~Ln- Hence we may write 
Tr [VU<t>)T N (ri\ = J2 (^V n N (cj>)T N (»)5 u ) 



umodN 



modiV V / V / 

m^H A/" \ / 



umodN 



for A" large and ||/z|| < A^ 1 / 2-5 , 5 > 0. As we shall justify below any smooth observable 
can be well approximated with a suitable cut-off in /x. In the semiclassical limit the matrix 
elements of U^(0) are given by fl8| 

(5 U , -^fY D (jj + m ^ ! Pl ex P [ 2 ^V + m l5 ^)] , (B.7) 

n mi=0 ^ ' L ^ ' - 

where 



dq 2 dqi 

S(f,n (q 2 , <7i) is the classical action of the lift of 4> n on the covering plane [10[|, and M n depends 
on the unperturbed map and on n. 

Inserting (|B.7|) into ( B.6|) and applying the Poisson summation formula yields 



n mi=0m 2 eZ J e 

x exp ^mAT (s^ (q + mi,q— - m 2 g) exp (27ri fx 2 q) dq (B.9) 
In leading order (N — > oo) we have, uniformly for all [i with < A^ 1 / 2-5 , 

£> (g + mi, g - ^) ~ D(g + m u q) (B.10) 



and 



a ( , a , , \ dS^(q + mu q' 
<V (g + mi , g -_) ~^"(g + mi,g) — 



<V(g + r%,g) +p(g)^ 



i'=i N (B.ll) 



Inserting ( p.lOj) into ( B.9|) gives 

M n »i_ e 



Tr [U^)T*(/*)] ^ivf ' exp [2th (/^(g) + /i 2 g)] x 



mi=0 rrt2£Z 

x D(g + mi, g) exp [27riAT {S(^{q + mi, g) - m 2 g)] dg. (B.12) 
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The above expression is the same as the formula that we would have obtained in de- 
termining the trace of the n-th propagator of an Anosov map with D n (qi,q 2 ) replaced 
by D n (qi, q 2 ) exp(27ri/Lt ■ z). Following the same steps which lead to the trace formula of 
Anosov maps |1| |10| we obtain 

Tr [UJ(*)T„(/»)] ~ V eXP( !^- Z ^ exp [2,riiV (S r (q, + m u q t ) - m 2 q,)} . (B.13) 

Let Jn be a truncated Fourier approximation to /, 

f N (z)= //xexp(27ri M -z). (B.14) 

|| M ||<ArV2-* 

Taking finite linear combinations, (|B.13|) yields 

N f ( ) 

YV-, Op N (f N ) Vj ) ~ V Ip^L exp [ 27 riiV (V(9f + m:, ®) - m 2 g f )] . (B.15) 



To conclude the proof of our claim ( |B.2| ), we note that, since the Fourier coefficients of / 
are decreasing faster than any power, both 

N N 

| YjiVj, Ov N (f)<Pi) ^ - ((Pi, Pjv (/ivM) e^l < £ Op N (f - f N )<pj)\ (B.16) 
and 

I E /(Zf) ,~ Mzt) ^ ^ N ( vfo + ^> ») - ^ gf )] | < E I /(Zf) /" M * f) I ( R17 ) 

are rapidly decreasing, as N becomes large. 
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